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Abstract 

We prove the global asymptotic stability of a well-known delayed negative- 
feedback model of testosterone dynamics, which has been proposed as a model 
of oscillatory behavior. We establish stability (and hence the impossibility of 
oscillations) even in the presence of delays of arbitrary length. 

Keywords: testosterone dynamics, monotone systems, negative feedback, global 
stability. 



1 Introduction 

The concentration of testosterone in the blood of a healthy human male is known 
to oscillate periodically every few hours, in response to similar oscillations in the 
concentrations of the luteinising hormone (LH) secreted by the pituitary gland, and 
the luteinising hormone releasing hormone (LHRH), normally secreted by the hy- 
pothalamus (see [n],|lH|)- In his influential textbook Mathematical Biology jlTj, J.D. 
Murray presents this process as an example of a biological oscillator, and proposes a 
model to describe it (pp. 244-253 in this edition). To obtain oscillations in an oth- 
erwise stable model, he introduces a delay in one of the variables, and by linearizing 
around the unique equilibrium point, he presents an argument to find conditions for 
the existence of such oscillations. This section in his book has remained virtually 
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unchanged since the first edition of 1989, up to the recent pubhcation of the third 
edition in 2002. 

The study of delayed models is one of great interest for its applicability in biolog- 
ical applications since it introduces a very relevant realism. (Consider for instance 
the delay between the moment a protein is transcribed, and the moment the folded 
and translated protein gets to act as a transcription factor back in the nucleus.) This 
realism often comes at the expense of a higher difficulty in mathematical treatment. 

As study" for a method for proving stability in a class of dynamical sys- 

tems with delays, we show in this paper that Murray's model in fact does not exhibit 
oscillations. The biological model itself, while simplified, is still interesting in its 
own right, and belongs to a commonly recurring class of models of negative feedback 
proposed (in undelayed form) by Goodwin 0, and illustrated in Goldbeter jBj. In 
what follows, we first study the linearized system around the unique equilibrium, 
estabhshing local stabihty, and then proceed to show the global stabihty of the sys- 
tem, borrowing ideas from monotone systems and the theory of control. We also 
propose an explanation for the confusion in fTI] . 

2 The Model, And Its Linearization 

The presence of LHRH in the blood is assumed in this simple model to induce 
the secretion of LH, which induces testosterone to be secreted in the testes. The 
testosterone in turn causes a negative feedback effect on the secretion of LHRH. 
Denoting LHRH, LH, and testosterone by R, L, and T respectively, and assuming 
first order degradation and a delay r in the response of the testes to changes in LH, 
we arrive to the dynamical system 



Here, bi,b2,b3, gi, g2 are positive constants, r > and f{x) = A/{K + x), although 
other positive, monotone decreasing functions could be employed as well (see Mur- 
ray, p. 246). 

By setting the left hand sides equal to zero, it is straightforward to show that 
there are as many equilibrium points of (Q) as there are solutions of 



R = f{T) - b,R 



L = giR — 62-^ 

f = g2L(t - r) - 63T 



(1) 




(2) 



namely for each such solution Tq of ((21), one has the equilibrium 




(3) 
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and by the assumption of positivity and monotonicity of / there always exists a 
unique solution of (j21), thus a unique equilibrium point of ((T)). Linearizing around 
that point we obtain the system 

X = f'{To)z - bix 

y = g^x- b2y (4) 
z = 92yit -r)- hz 

The characteristic polynomial of which determines all solutions of (0} of the 
form v(t) = VqC'^*, is 

{\ + h){\ + h){\ + h) + de-^^ = d=-f\To)gig2>0 (5) 

Proposition 1 The linear system ^ is stable, for all values of bi,b2,b3, gi, g2,T 
and f{x) = A/{K + x). 

Proof: For there to be an unstable solution of (jlj), there must be a solution A of 
(jSj) such that Re A > 0. Assuming that this is the case, we have 

d> \-de-^^\ = IA + 61IIA + 62IIA + 63I > IhlMM =6i62&3- (6) 

But on the other hand, using the choice for /(T) above, we have f'{To) = —Aj (K ^ 
To)' = -fiTo)/{K + T),and 

d = -f(To)gig2 = IJ^^^l, 9i92 = &i&2&3^7^^ < ^i&2&3 , (7) 
which is a contradiction. Q.E.D 



3 Global Asymptotic Stability of the Model 

Even with the addition of only one simple delay, it is probably best to view (^J) as 
a dynamical system with states in the space X of continuous functions from [— r, 0] 
into the closed positive quadrant M.^. The right hand side of (P) defines a function 
F : X ^ M.^ in the natural way, and given an initial state G X, the solution 
of the system is the unique absolutely continuous function x : [— r, 00) such 
that 

X(o) = and x(t) = F(x(t)), t >0 (8) 

Here, X(t), or simply Xt, is the state 7(5) = x(t + s), s G [— r, 0]. The function 
0) = Xt will be from now on formally identified with system For proofs of 
the fact that $ is well-defined, and more details, the reader is referred to |^ ISllTlj. 
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Cutting the Loop We define a funcion G : X x M.^ in a very similar 

manner to F: for 0(s) = (-R(s), L(s), T(s)), let 



Given a piecewise continuous function u : M+ — > M+, called an input , we define 
m) = Xt, where x : [— r, oo) M.^ is the unique absolutely continuous func- 
tion such that 



In effect, we are thus cutting the feedback loop induced by T upon R, and replacing 
it with an arbitrary input u{t). 

Notation: given x,y E R'^, let x < y denote Xi < yi, i = 1,2,3. For (f),ip E X, 
let (f) <iIj denote <^p{s), Vs G [-r, 0]. 

Theorem 1 The dynamical system with input "^{t, (p, u) satisfies the following prop- 
erties: 

1. If the input u{t) converges to w E M+, then '${t,(j),u) converges as t tends to 
oo towards the constant state 



for any initial state G X. 

2. Let ui,U2 be inputs, and pick any two initial states 0, ■?/' G X. // Ui{t) < 
U2{t) Wt and (j) < ip, then 0, Ui) < \E'(t, ip, U2) Vt. 

Proof: Suppose that u{t) converges towards w G M+, and let G X arbitrary. 
The dynamics of the component R{t) of the solution x{t) is determined by the 
equation R{t) = u{t) — biR{t), and so R(t) converges towards w/bi. Applying a 
very similar argument to L{t) and T{t) in this order, we obtain the first result. 

The proof of the second statement follows by the "Kamke condition" (see |14j): 
if wi < W2, (p < ip, and 0(O)i = ip{0)i (that is, the ith components of and ip 
are equal), then G{(p,wi)i < G{ip,W2)i. For instance, if = (_Ri,Li,Ti), ip = 
{R2,L2,T2), < V, and i?i(0) = i?2(0), then wi - 6ii?i(0) < W2 - 6ii?2(0). This 
can be checked for L and T in the same way. The fact that the Kamke condition 
implies the desired property follows from the results in jHj; however, in the interest 
of exposition and since the proof is so short, we provide it next. 

-'^We won't require the more general control-theoretic definition where u is measurable and 
locally bounded, see P 



F{(P,w) = {w- b,R{0), g,R{0) - 62^^(0), g2L{t - r) - bsT{0)) . 



X(o) = and x{t) = G{xt, u{t)), t > 0. 



(9) 
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Let x(t) be the solution of Q with input Ui and initial condition 0, and let 
= G + (e, e, e), for e > 0. Let ye(^) be the solution of y(t) = G^iyt, U2) with initial 
condition ifj. Suppose by contradiction that at some point ti, x{ti) ^ Veiti), and so 
there exists a component i (that is, R, L or T) and to such that xt^ < x(to)i = 
ye{to)i and x(to)i > ye{to)i- But then 

i;()f:o)i = G{xto,Ui{to))i < G{y,to,U2ito))i < G^{y,to,U2{to))i = y,(to)i 

which is a contradiction. We thus conclude that x{t) < y^{t), Wt > 0. Now, it 
can be shown (jH|,Pl|) that as e ^ y^{t) converges pointwise to yit), the solution 
of Q with input U2 and initial condition ip, and from here the conclusion follows. 
Q.E.D. 



Definition 1 Given x : [— r, 00) M.^ be an arbitrary trajectory, we say that 
z G M.^ is a lower hyperbound of x{t) if there is zi, Z2, ■ ■ ■ — ^ z and ti < t2 < 
t^ ...—>■ 00 such that for all t > ti, Zi < x{t). A similar definition is given if for all 
t >ti, Zi> x{t), and we say that z is an upper hyperbound of x{t). 

For instance, z is a lower hyperbound of the trajectory x if it bounds from below 
x{t) for every t. Similar definitions are given for inputs u{t). The previous Theorem 
is the basis for the following result. 

Theorem 2 Let v G M+ be a lower hyperbound of the input u{t), and let cf) E X be 
arbitrary. Then k{y) is a lower hyperbound of the solution x{t) of the system 
If V is, instead, an upper hyperbound of u{t), then k{v) is an upper hyperbound of 
x{t). 



Proof: Suppose that f is a lower hyperbound of u{t), the other case being similar, 
and let f 1, f2, . . . — > f and ti < t2 < • • • ^ C)0 be as above. 

For every z > 1, let ?/j G Mj^ and Vi C M?^ neighborhood of k{vi) that is open in 
M"^, chosen in such a way that yi<Vi and \yi — k{vi)\ < Also, let 



Ui{t) 



u{t), 0<t <tr. 



Let Ti < T2 < ... 00 be defined by induction as follows: Ti = 0, and if Tj_i is 
defined, let Tj be chosen such that Tj > Ti_i, Ti > ti, and for all t > Ti : Xi(t) = 
\E'(t, 0, Mj) G Vi. By the previous theorem, Xi(t) < x(t)Wt, and so yi < x(t), Vt > Tj. 
As yi — > k{v), the conclusion follows. Q.E.D. 

The following simple Lemma is standard in the literature on discrete iterations 
(and is used in a similar context in 3J); we provide a proof for expository purposes. 
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Lemma 1 Let S : M+ be a continuous, nonincreasing function. Then the 

discrete system Un+i = S{un) has a unique, globally attractive equilibrium if and 
only if the equation S{S{x)) = x has a unique solution. 



Proof: If the system has a unique, globally attractive equilibrium u, then this 
point is a solution of the equation S'^(x) = S{S{x)) = x. Any other point u cannot 
be a solution of this equation, as S"{u) must converge to u. This proves one of the 
directions of the lemma. 

Conversely, suppose that the equation S'^(x) = x has a unique solution. Let 
u G M+ be arbitrary, and consider the sequence = S"'{u). If u < U2, then since 
S"^ is a nondecreasing function, we have U2 < u^, and so 

U < U2 < < Uq < . . . . 

But the sequence M2,'U4... is bounded (by 5'(0)), and so U2n must converge to 
some point vq. The same argument applies if U2 < u, and also for the sequence 
Ml, U3, M5, . . ., which must converge to some point f 1. But the continuity of S implies 
that both vq and Vi are solutions of S^{x) = x, so Vq = Vi are both equal to our 
unique solution, and thus converges to this point, independently of the choice of 
u. Q.E.D. 

Consider for instance S{x) = p/ {q + x), where p, q are positive real numbers. If 
X satisfies S'^(x) = x, then it holds that 

X = ^ 
q + S{x) 

which can be rearranged as x'^ + qx—p = 0. Using the quadratic formula, it becomes 
clear that there is always exactly one positive solution. 
This example will be useful in what follows. 

Theorem 3 All solutions of the system with f = A/{K + x), converge towards 
the unique equilibrium, for any choice of the parameters bi,b2,b3, gi, g2,T, A, K . 



Proof: Consider any initial condition G X, and the corresponding solution 
x{t) = {R{t),L{t),T{t)) of dH Defining the input u{t) = /(T(t)), and using it to 
solve the system with initial condition 0, we arrive of course at exactly the same 
solution x(t). 

Let V bound u{t) from below for all t - for instance, v = will do. Then by 
Theorem 121 k{v) is a lower hyperbound of x{t). In particular. 
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is a lower hyperbound of T{t). But, since / is a nonincreasing function, this implies 
that f{Qv) is an upper hyperbound of f(T(t)) = u(t). Defining Vi = f{Qv), we 
apply the same theorem once again to show that k{vi) is an upper hyperbound of 
x{t), t>2 = fiQvi) is a lower hyperbound of u{t), etc. But 

K + Qx q + X 

for p = A/Q, q = K/Q. Thus we see that f„ = S"'{v) is a convergent sequence of 
numbers that are alternatively upper and lower hyperbounds of u{t). This easily 
implies that u(t) itself is a trajectory that converges to the unique solution u of the 
equation S'^{x) = x. By Theorem^ x{t) converges towards k{u), independently of 
the choice of the initial condition 0. 

Finally, this implies that k{u) is the unique equilibrium of the system, other- 
wise one could reach a contradiction by taking this equilibrium as constant initial 
condition. Q.E.D. 



3.1 Discussion 

Several remarks are in order: first, the value of the delay r was almost never used, 
and indeed can be arbitrarily large or small. In fact, we can introduce different 
delays, large or small, in all of the first summands of the right hand sides of (jH}, 
and the results will apply with almost no variation. If delays are introduced in 
the second summands, the system will not be monotone, that is, won't satisfy the 
second property of Theorem ^ which is essential for this argument. But then 
again, introducing a delay in the degradation terms wouldn't be very biologically 
meaningful. For more on monotone systems, the reader is referred to the excellent 
textbook by Hal Smith [T3], and jU]. 

The above argument is an illustration of a more general treatment on a class 
of delayed dynamical systems with monotone subsystems and negative feedback in- 
terconnection. The underlying order may be generalized as x < y iS y — x lies in 
a cone K (1 M.^ (see ^Q). This provides for more generality and applicability in 
biological problems. The key sufficient condition is that the discrete dynamical sys- 
tem Un+i = S{un) be globally attractive; in a sense the dynamics of the continuous 
system is reduced to that of the discrete one, which may eventually involve state 
spaces with substantially fewer dimensions. See PP,|21 and work to appear by the 
present authors for this more general treatment. 

As for the conclusions in pp. 244-253 of Mathematical Biology, we may venture 
to suggest that in eq 7.49, p. 247, the author writes the characteristic equation (jSj) 
of the linearized system (0)) as 

\^ + a\'^ + b\ + c + de-^^ = (10) 
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where a,b,c,d are all written in terms of the original parameters of the system: 
a = 61 + 62 + ^3) etc. From here on the efforts are concentrated in finding a root A 
of this equation with Re A = 0, for some well-chosen coefficients a, b, c, d. But the 
author seems to disregard in the remaining argument the fact that a, 6, c, d cannot 
be chosen arbitrarily and independently, hut rather that their values are determined 
from choosing arbitrarily &i, &2; ^s? fi'i? fi'2; t- Thus for instance, it is assumed in the 
last line of p. 251 that d > c, without justification from the original variables. The 
former assumption turns out not to be possible to satisfy for the particular choice 
of /, as seen in the proof of Proposition ^ 

We point out that a simple modification can make oscillatory behavior possible. 
In p. 246 of ^T], the author discusses varying cooperativity coefficients of f{x) = 
A/{K + x"^), then settles for m = 1 for the delayed model. If indeed m is increased, 
then it is very possible to have d > c and the remaining argument in the section 
will be valid. One example of this is when parameters are picked as follows: 



Another interesting contribution to the modeling of testosterone dynamics is the 
paper [T2j by Ruan et al., where sufficient conditions are found for stable and oscil- 
latory behavior in a neighborhood of an equilibrium. We would like to describe the 
relationship between fT^J and our own result, given the similarity of the hypotheses 
and the potentially conflicting conclusions: global stability in our results vs. Hopf 
bifurcations in ^2]- Moreover, we will simplify the statement of that result. In 
that paper, several new quantities are introduced in order to state the main result. 
Theorem 3.1. In terms of the original variables of the system (61, 62, &3, etc.), these 
are as follows: 



m = 2, A = W, K = 2, b 



'I = 



1, &2 



1, (71 = 10, 92 = 10. 




Theorem 3.1 holds under the assumption that 



(fel + &2)(&l+fe3)(fe3 + &2) <d 



(11) 



and deals essentially with three following three special cases: 



1. 6162&3 > and A < 0, 



2. 6162&3 > d and zi > 0, 
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3. < d, . 

In case 1, (local) asymptotic stability is guaranteed for arbitrary delay lengths (part 
(i) of the Theorem), while in cases 2 and 3, and under some additional conditions 
(parts (ii) and (iii) of the Theorem), stability holds for small enough delays, but a 
Hopf bifurcation occurs at some critical value of this delay length. In light of the 
above computation, case 1 can never be satisfied (for variables p, q, r generated 
from the original set of parameters bi, 62, etc.). Similarly, the condition zi > 
will never be satisfied, since 

zi > <^ A > <^ 3q < 

so case 2 cannot hold either. One is only left with case 3, which is actually a 
consequence of (II Ij) . On the other hand, for the particular choice of f{x) made in 
[TT] and the present paper. Proposition ^ shows that we always have 6162^3 > d- 
Thus Theorem 3.1 does not apply for the present model, as well as for any choice 
of the function / and any set of parameters such that 6i62^'3 > d. 

Acknowledgement: We would like to thank Augusto Ponce for useful sugges- 
tions. 
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